Salt-induced collapse and reexpansion of highly charged flexible polyelectrolytes 
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We study the salt-dependent conformations of dilute flexible polyelectrolytes in solution via com- 
puter simulations. Low concentrations of multivalent salt induce the known conformational collapse 
of individual polyelectrolyte chains, but as the salt concentration is increased further this is followed 
by a reexpansion. We explicitly demonstrate that multivalent counterions can overcompensate the 
bare charge of the chain in the reexpansion regime. Both the degree of reexpansion and the occur- 
rence of overcharging sensitively depend on ion size. Our findings are relevant for a wide range of 
salt-induced complexation phenomena. 
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Many biological and synthetic polyelectrolytes undergo 
two macroscopic phase transitions upon addition of mul- 
tivalent salt or charged small molecules 0, 0, E3 [3 IK 
First, precipitation of polyelectrolytes occurs if the con- 
centration of the added salt exceeds a critical value. This 
precipitate subsequently redissolves if the concentration 
is increased beyond a second critical value. These phe- 
nomena, jointly referred to as reentrant condensation 
have attracted considerable attention because they are 
a fundamental and generic aspect of polyelectrolyte be- 
havior, with potential relevance for the understanding 
and development of biological phenomena and applica- 
tions such as gene delivery Q. Under dilute conditions, 
the formation and redissolution of multimolecular ag- 
gregates are replaced by, respectively, single-chain "col- 
lapse" and reexpansion. This behavior is directly ob- 
served in, e.g., the stretching of individual DNA coils in 
salty solution 0. 

Whereas several theoretical explanations for reentrant 
condensation have been proposed 0, 0, Hoi HlT |. im- 
portant open questions remain regarding the suggested 
mechanisms, and several predictions have only partially 
been verified. Two main scenarios can be distinguished. 
(1) Charge inversion: Counterions form a strongly corre- 
lated liquid at the polyelectrolyte surface and, at high salt 
concentration, overcompensate the bare chain charge. 
Condensation of a polyelectrolyte is predicted to take 
place if its effective charge is close to zero [(jllSj. (2) The 
two-state model |I| explains precipitation and redissolu- 
tion of flexible polyelectrolytes by assuming that individ- 
ual chains adopt either a collapsed or an extended state. 
This scenario predicts a strong sensitivity to ion size, 
which is not considered in Refs. 0,01- I n addition, in 
this case overcharging does not occur by necessity, but 
may be mitigated by coion association 

Computer simulations offer the possibility to study this 
behavior at a microscopic level, but hitherto most simu- 
lations of flexible chains have focused on solutions with 
only divalent salt [3 Q or no salt at all [H HEl E3. 
Whereas chain collapse has been observed under the in- 



fluence of divalent salt or counterions |16| . the reexpan- 
sion observed in experiments and predicted by theory has 
not been reproduced in computer simulations. In this 
Letter, we aim to clarify this apparent discrepancy and 
to understand the behavior of strong polyelectrolytes in 
solutions of multivalent salt by studying the conforma- 
tions of dilute chains over a wide range of salt concen- 
trations. Our calculations allow us to explicitly address 
the theoretical predictions. Specifically, we answer the 
following questions: (i) How does the conformation of a 
polyelectrolyte depend on salt concentration and valency, 
and under what conditions can an extended (i.e., non- 
collapsed) structure be observed at high concentrations? 
(ii) Is there a dependence on ion size, as suggested by 
some theories? (iii) Can the charge of a polyelectrolyte 
be overcompensated by condensed ions, such that its ef- 
fective charge is reversed compared to the bare charge? 

Following Stevens and Kremer jl5j, we model an an- 
ionic polyelectrolyte as a bead-spring model in a continu- 
ous medium of uniform dielectric constant e, representing 
the solvent. A chain consists of N monomers of size a, 
which each carry a charge — e. Adjacent monomers are 
connected by a FENE bond potential 15] with maximum 
extension Rq = 2a and spring constant k = 7£lj/<t 2 . 
The counterions are monovalent, whereas the added salt 
dissociates into Z-valent cations and monovalent anions. 
Electrostatic interactions are computed via the Ewald 
summation and the excluded volume of monomers and 
ions is modeled via the shifted-truncated Lennard- Jones 
potential with coupling constant £lj = k^T/1.2, where 
&b denotes Boltzmann's constant and T the absolute 
temperature. Periodic boundaries are imposed. We ex- 
press the Coulomb interaction between particles i and j 
as k^TX^ZiZj Jt, where Ab = e 2 / (Airseok-BT) is the Bjcr- 
rum length, Eq the vacuum permittivity, and Zi and Zj the 
valencies of the particles. In aqueous solution at room 
temperature, Ab = 7.14 A. We set Ab = 3ct, represen- 
tative of a prototypical flexible polyelectrolyte such as 
sodium poly(styrene sulfonate). 

To compute conformational properties, we employ 
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FIG. 1: Squared radius of gyration Rg of a polyelectrolyte of 
length TV = 64 as a function of salt concentration C s . Curves 
corresponds to added salt with counterion valency from 1 to 4. 
The curve 3:1* refers to a system with ions that have a twice 
smaller diameter (see text for details). The dashed line cor- 
responds to a reference system without any electrostatic in- 
teractions. Error bars are smaller than the symbol size. 



Monte Carlo simulations, which makes it possible to 
accelerate the relaxation of longer chains via repta- 
tion moves. Since we focus on the dilute polyelec- 
trolyte regime, relatively large simulation cells are re- 
quired, resulting in considerable ion numbers in sys- 
tems with added salt. For computational efficiency, we 
therefore confine ourselves to a monomer concentration 
C m = 0.008cr~ 3 , which corresponds to 0.99mol/£. This 
is higher than typical experimental conditions Q, but 
we have explicitly verified that our findings are repre- 
sentative for the entire concentration range 0.001cr~ 3 < 
C m < 0.008er~ 3 . In addition, for all chain lengths stud- 
ied here (N = 16, 32, 64, 96), C m lies well below the 
overlap threshold C m as estimated from the radius of gy- 
ration R g . We adopt cr, <r~ 3 , and e as our units of length, 
concentration, and charge, respectively. 

We first study the conformation of a polyelectrolyte of 
length N = 64 in the presence of (Z:l)-salt, where Z is 
varied from 1 to 4. Ions and monomers have the same 
size. Figure Q] shows the average squared radius of gyra- 
tion Rg as a function of salt concentration. In the absence 
of salt, the chain adopts an extended conformation, owing 
to the electrostatic repulsions between monomers. Upon 
addition of monovalent salt, these repulsions are screened 
and Rg gradually decreases. By contrast, for multivalent 
salt, a much stronger decrease in Rg is seen, occurring at 
considerably lower salt concentrations. This is the con- 
formational collapse that has previously been observed 
for polyelectrolytes with multivalent counterions under 
salt- free conditions 0, 0|. The smallest value of R 2 g 
occurs for C s near Cz, the (Z:l)-salt concentration at 
which the total charge of the Z-valent cations neutralizes 



the bare polyelectrolyte charge. Accordingly, this com- 
pact state occurs at a salt concentration that decreases 
with increasing valency Z , consistent with the two-state 
model 0, Il0| . In addition, higher counterion valency 
leads to more compact states. A striking effect occurs 
once the salt concentration is increased beyond Cz- The 
chain starts to swell again, in close analogy with the re- 
dissolution observed for multichain aggregates 

0,13 Al- 
though the chains that are collapsed by trivalent and 
tetravalent counterions reexpand strongly, they remain 
more compact than the chain in 2:1 salt for all salt con- 
centrations investigated. Comparison with the structure 
of the same polyelectrolyte in monovalent salt, which ex- 
hibits a slow, monotonic decrease of Rg with increas- 
ing C s , emphasizes the important role of counterion va- 
lency. The salt concentration at which the reexpanded 
state is reached is insensitive to monomer concentration, 
whereas the collapse concentration Cz is proportional 
to C m . Thus, our relatively high choice of C m causes 
both conformational changes to occur over a narrower 
range in salt concentration than in typical experiments. 

We further quantify the swelling of the polyelectrolyte 
coil by investigating the variation of R g with degree 
of polymerization N, for tetravalent salt. As antici- 
pated from scaling theory 19], the radius of gyration 
is described by a power law N v for the range of chain 
lengths studied here (16 < N < 96). Figure de- 
picts the concentration-dependent exponent v extracted 
from least-squares fits at fixed C s . At low concen- 
trations, v exceeds the self-avoiding random-walk value 
^saw = 0.588, owing to the electrostatic repulsions be- 




FIG. 2: Scaling dimension v of the radius of gyration for 
a polyelectrolyte in a tetravalent salt solution. At low salt 
concentration C s , v exceeds the self-avoiding random- walk 
value 0.588, owing to the electrostatic repulsions between 
monomers. For the concentration Cz=4 corresponding to the 
minimum in Fig.0 the exponent reaches the value for a com- 
pact globule. The subsequent reexpansion of the chain at 
higher salt concentrations is reflected in an increase in v. 
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FIG. 3: Net charge Q(r) within a tube of radius r around a 
chain (N = 64) at various salt concentrations C s of trivalent 
salt, (a) For ion diameter a (case 3:1 in Fig.0 an overcharg- 
ing peak appears at short r. (b) For smaller ions (case 3:1*) 
the overcharging peak is absent. 



tween monomers |20| . Near the neutralization concentra- 
tion Cz=4 = 0.002, the chain attains its most compact 
state and v — 0.34 ± 0.02, in good agreement with the 
poor-solvent value, v = i. Upon reexpansion, the expo- 
nent v increases and approximately reaches vsaw at the 
highest concentrations. 

Chain reexpansion has been linked to charge reversal 
(or overcharging) [l2^ . Since for spherical colloids over- 
chargingwas found to strongly depend on ion excluded 
volume |2l|, we proceed to clarify the effect of ion size 
in the case of flexible polyelectrolytes. To this end, we 
repeat our calculations for 3:1 salt in which all ions have 
a diameter a/2, i.e., twice smaller than the monomers. 
The electrostatic attraction between a trivalent counte- 
rfoil and a monomer at contact thus increases by a factor 
| and equals the attraction experienced by tetravalent 
ions of diameter a. As the smaller ions allow the chain 
to adopt a more compact structure, it may appear plau- 
sible that in this situation R g is even smaller than for 
tetravalent salt and hence also smaller than R g for triva- 
lent salt of diameter a. However, as shown in Fig. ^ 
(cf. the curve labeled 3:1*), the actual situation is more 
complicated. The polyelectrolyte chain is indeed less ex- 
tended than in tetravalent salt solution for low salt con- 
centrations and shrinks further if C s is increased. But 
this decrease is weaker than for larger ions, and near the 
neutralization condition the chain is much less compact 
than for tetravalent or even trivalent salt of diameter a. 
For higher salt concentrations we observe a gradual ex- 
pansion of the chain, and R g reaches a magnitude that 
is approximately the same for both trivalent cases. 

The mechanism for this remarkable dependence on 
ion size can be understood from the charge distribu- 
tion around the polyelectrolyte. We define a worm-like 
tube around the chain, consisting of N spheres of ra- 



dius T centered around the monomers, and monitor the 
ensemble-averaged net charge Q(r) within this tube, as a 
function of salt concentration. As shown in Fig. E^a) for 
trivalent salt with ion size a, Q(r) is always negative for 
< C s < Cz=3 = 0.00267 and only approaches zero at 
large r, to satisfy electroneutrality. In this regime, upon 
addition of salt all trivalent counterions condense on the 
chain and replace the monovalent counterions, as pre- 
dicted || . As C s approaches Cz=3 the magnitude of the 
net charge Q{r) decreases, until the chain is completely 
neutralized by the condensed trivalent counterions and 
takes a compact structure. For C s > Cz=3 , not all triva- 
lent ions condense, although the number of condensed 
trivalent counterions can exceed the number required to 
neutralize the chain, in accordance with Refs. 0, fl2j |. 
This overcharging is reflected in the appearance of a 
strong peak at small tube diameter r in Fig. [SJa). The 
excess charge attracts coions that form a second layer, 
resulting in a rapid decrease of the peak with increas- 
ing r. At the highest salt concentration investigated, 
Q(r) exhibits an oscillatory tail, indicative of a multi- 
layered ionic structure, similar to what has been observed 
for electrical double layers , colloids |2^| and rod-like 
polyelectrolytes Interestingly, for monovalent salt 

(not shown) no overcharging peak occurs, consistent with 
the absence of important conformational changes when 
the salt concentration is varied. 

The behavior of Q(r) drastically changes upon reduc- 
tion of the ion size, as illustrated in Fig. |3Jb) for case 
3:1* (ion diameter cr/2). Again, the net charge within the 
tube decreases with increasing salt concentration, but it 
is not neutralized at Cz=z, consistent with a much larger 
coil size than for the system with larger ions (Fig.QJ. For 
C s > C'z=3 a small peak appears in Q(r), but there is 
hardly any overcharging even at the highest concentra- 
tions. A study of the ion configurations in the simula- 
tions reveals that this remarkable difference arises from 
the strong increase in counterion-coion attraction if their 
size is reduced. This association, akin to the formation of 
Bjerrum pairs, reduces the effective valency of the salt. 
Consequently, compared to the 3:1 case, the chain ex- 
hibits a less strong conformational collapse and Q(r) re- 
mains more negative. When pairs of trivalent counterions 
and monovalent coions condense on the chain, they tend 
to be oriented with the trivalent ion toward the monomer, 
explaining the appearance of a small peak in Q(r). Anal- 
ogous behavior has been observed for simple colloids [2l| , 
indicating that the decrease in overcharging upon reduc- 
tion of ion size is a generic effect. It also explains why 
neither reexpansion nor overcharging of a polyelectrolyte 
were observed in earlier simulations of a flexible chain in 
the presence of divalent salt ^ij. In those simulations 
the ion sizes were chosen as half the monomer size |2^ . 

It is tempting to assign an "effective" charge to the 
chain based upon the net tube charge Q(r). However, 
as Fig. |2| indicates, such an interpretation is arbitrary if 
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FIG. 4: Drift velocity Vdrift ( m units of cr/r) of the center 
of mass of a polyelectrolyte (N = 16) as a function of salt 
concentration C s , in the direction of the electric field. The 
sign reversal in the electrophoretic mobility provides direct 
evidence for charge reversal by multivalent counterions. 



oscillatory behavior in Fig- EI 27]. 

In summary, we have demonstrated, for the first time, 
that dilute flexible polyelectrolytes in a solution with 
multivalent salt undergo not only a conformational col- 
lapse, but also a reexpansion. These findings support the 
two-state model for polyelectrolyte conformations 0, 0] 
and confirm the important role of ion excluded vol- 
ume Furthermore, we have presented unambiguous 
evidence for overcharging by means of condensed multi- 
valent counterions at high salt concentrations 0, 0] , in 
which the reversal of the effective charge coincides with 
the most compact chain conformation. 

This material is based upon work supported by the 
U.S. Department of Energy, Division of Materials Sci- 
ences under Grant No. DEFG02-91ER45439, through 
the Frederick Seitz Materials Research Laboratory at the 
University of Illinois at Urbana-Champaign. 



there is no clear prescription for the proper choice of the 
tube radius. To circumvent this ambiguity, we study the 
electrophoretic mobility as a function of salt valency and 
concentration for a polyelectrolyte of TV = 16 monomers 
via Langevin dynamics simulations, using the LAMMPS 
package [2^. The collision frequency in the thermostat 
equals 7 = r _1 , where r = a^J m/(k^T) is the time unit 
and m the ion mass. We apply a constant electric field of 
magnitude E — Q.lk^T / (ecr), which, as we carefully veri- 
fied, does not distort the conformational pro per ties of the 
chain or the charge distribution around it [27J. Due to 
the absence of hydrodynamic effects this calculation does 
not realistically reproduce dynamic properties, but it will 
properly reflect the sign of the electrophoretic mobility 
and can thus be used to probe the effective charge with- 
out imposing an arbitrary criterion for the condensation 
of specific ions. Indeed, the center-of-mass drift velocity 
of the chain (Fig.0} reveals a striking dependence on salt 
concentration. For both trivalent and tetravalent salt, it 
exhibits a zero crossing, i.e., the electrophoretic mobility 
of the chain changes sign upon addition of salt, as ob- 
served experimentally for colloids 28] . This sign reversal 
occurs close to the neutralization condition (minimum in 
Fig. IT]) and thus depends on counterion valency. Further- 
more, for trivalent salt with smaller ion size (case 3:1*) no 
sign reversal occurs, i.e., no overcharging takes place. We 
view these observations as a strong confirmation of the 
charge reversal suggested by Fig. and as evidence for 
the relation between charge reversal and reexpansion, but 
also note that a less strong collapse and subsequent reex- 
pansion can occur without overcharging. At the highest 
salt concentrations, a second sign reversal occurs, which 
is possibly related to the formation of a multi-layered 
structure of counterions and coions, consistent with the 
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